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We study the effect of cosmological expansion on orbits-galactic, planetary, or atomic-subject to 
an inverse-square force law. We obtain the laws of motion for gravitational or electrical interactions 
from general relativity-in particular, we find the gravitational field of a mass distribution in an 
expanding universe by applying perturbation theory to the Robertson- Walker metric. Cosmological 
expansion induces an (a/a)r force where a(t) is the cosmological scale factor. In a locally Newtonian 
framework, we show that the (a/a)f term represents the effect of a continuous distribution of 
cosmological material in Hubble flow, and that the total force on an object, due to the cosmological 
material plus the matter perturbation, can be represented as the negative gradient of a gravitational 
potential whose source is the material actually present. We also consider the effect on local dynamics 
of the cosmological constant. We calculate the perihelion precession of elliptical orbits due to 
' the cosmological constant induced force, and work out a generalized virial relation applicable to 

D gravitationally bound clusters. 

PACS numbers: 04.90.+e, 98.80.Jk 

(N ■ 

<n : 

7— I . 

i. introduction 

\q ; 

,_H ! If the universe has expanded by a factor of 1000 since the first atoms were formed, and by 30 some percent since the 
£NJ ■ formation of the solar system, it seems reasonable that there must be some effect, even if small, of this expansion on 
' the orbits of atoms and planets. On the other hand, spacetime is locally flat, and first-order deviations from flatness 
are accounted for by Newtonian gravitational theory, leaving no room for stray cosmological effects on local systems. 
So, do atomic or planetary orbits respond at least a bit to the universal expansion or don't they? Our purpose in 
this article is to answer this question both within matter only cosmological models and in the currently accepted 
O" 1 benchmark model including a cosmological-constant derived dark energy component. 

, Many authors have studied the effects of cosmological expansion on local dynamics with various and often conflicting 
conclusions. Most of the early work weighs against any cosmological effect on planetary orbits. A seminal contribution 
was made by McVittie Q, who found an exact solution to Einstein's equation that he interpreted as representing 
a point mass embedded in a Robertson- Walker expanding universe. McVittie, and later Jarnefelt @, Q found that 
planetary orbits in such a model do not participate in the cosmic expansion. Recent work on problems and progress 
in the interpretation of McVittie's solution is summarized by Nolan who considers the influence of the cosmic 
expansion on planetary orbits in the McVittie universe still to be an open question, fn another influential study, 
Einstein and Straus Q (see also Schiicking [(|) showed that a truncated Schwarzschild solution could be embedded 
in a background expanding universe with smoothly matched boundary conditions. In this picture, orbits around 
the central mass are completely unaffected by the expansion. However, this approach has been criticized on various 
counts, summarized by Carrera and Giulini (tJ. More recently, Dicke and Peebles argued on general grounds that 
there are no local effects Q . Their work was criticized by Noerdlinger and Petrosian Q and by Carrera and Giulini 0] 
on the grounds that they neglected the "friction" term in their equation that would have produced a positive effect. 
Anderson [To| obtained an equation of motion including such a friction term, but found that there are no expansion- 
induced effects on circular orbits, at least in the leading order of approximation. On the other side, many authors 
find small but non- vanishing effects on local orbits due to cosmological expansion. @, H, EH, EE EE EE EE EE E3 We 
will discuss many of these results in the course of this work. 

We base our approach on relativistic cosmology, and show that locally this reduces to a Newtonian framework. 
For local systems this is an appropriate approximation, as has been discussed by McCrea and Milne [IE [lj| and 
particularly by Dicke and Peebles Q and Callan, Dicke, and Peebles 20] for the present application. We restrict 
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our attention to a spatially flat universe, which is in accord with recent observations. In any case, a flat space 
approximation will be adequate for space scales much less than the radius of curvature. We consider models involving 
matter and dark energy only, and neglect radiation as of negligible import in the present epoch. 



II. EQUATIONS OF MOTION INCLUDING COSMOLOGICAL EFFECTS 

In this section we derive the equation for non-relativistic motion of a point particle in the electric or gravitational 
field of a spherically symmetric static source. The equation has the same form for either case: 

■± 01 „ a _ 

r = -r H — r (1) 

r z a 

where r is the position of the point particle in physical coordinates, a represents the strength of the interaction, and 
a(t) is the cosmological scale factor. 

The provenance of ([1]) is quite different for the two cases. For electric attraction, it comes from the curved space 
generalization of the Lorentz force law 

x»+r*x a xP = ±F» V V, (2) 

* 771 

where r£o are the Christoffel coefficients, F^ v is the electromagnetic field tensor, x^ is the particle coordinate in 
co-moving coordinates (so r — a(t)x), and the dot denotes a derivative with respect to proper time. On the other 
hand, for gravitational attraction, the equation of motion is the geodesic equation 

as" + K /3 x a x f3 = (3) 

with no separate force term-the effect of the source is contained in T^g. Since the electric case is more straightforward, 
we will consider it first. 

The Robertson- Walker metric for a flat expanding spacetime has the form 

^ = d -a\t)Si) ' (4) 



where a(t) satisfies the cosmological Einstein equations 
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= 8irGp . 
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~-2- = SnGp, (5) 
a z a 



in terms of the total cosmological proper density p and pressure p. The corresponding Christoffel coefficients 

r a/3 = (9ua,3 + 9ufi,a ~ 9 a 0,u) (6) 

are = aaSij, Tq = r* = (d/a)<5y, with all others vanishing. The field tensor is 



' ij ~ uuu « ' Oj ~ x jO ~ W 

r ( 

K —Ei —eijkBk 

where Et and Bi are related to the electric and magnetic fields. The field tensor satisfies Maxwell's equations 



F^. v = inj" , (8) 



For a spherically symmetric static source, Maxwell's equations require B% — and 



& = ^ (10) 
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outside the source, where x — \/ x\ + x\ + x\ , f — x/x, and C is some constant. Now the spatial components of the 
Lorentz force law (|2|) lead to 

x i + 2-x i x a = 2--^-f i x a . (11) 



For non-relativistic motion this is 



' -4ti'--- ( 12 ) 



a m a°x 

and in terms of the physical coordinates r = a(t)x it takes the form 

a Qqh 

n — n = it (13) 

a m r z 

since we can now identify C with the source charge Q. 

We derive the equation of motion for a test mass in a gravitational field in two steps. First, we use the Einstein 
equation to find the metric for an expanding universe containing a localized non-relativistic distribution of matter 
in addition to the uniformly distributed cosmological material. We restrict ourselves to small distortions from the 
cosmological background. After obtaining the metric we use the geodesic equation in the weak-field limit to obtain 
the equation of motion. 

A massive object such as the sun inserted into a background expanding universe perturbs the usual Robertson- 
Walker metric. Exact solutions describing inhomogeneous cosmologies have been extensively studied (reviewed for 
example in the book of Krasihski pUp. but for our purposes a perturbative approach will be adequate physically 
and more straightforward to interpret. Accordingly, we will consider perturbations to the Robertson- Walker metric. 
[22l [23L l24l . [2{| For the weak-field case with only non-relativistic motion of localized source material it is sufficient to 
consider the scalar perturbation 

9 ^-{ -(l-2 X (t,x))a 2 (t)S t] ) [M) 

and a stress-energy tensor of the form 

T » v = X! + VaW^Uv - VaQnv] + PU^U U (15) 

a 

where p a and p a are the proper density and pressure of the a th component of cosmological material (matter: m, dark 
energy: A) and p(t, x ) is the proper density of the matter perturbation. To emphasize, p m + p represents all of the 
matter actually present, so p itself is the actual (total) matter minus the cosmological matter. The four- velocity 
has only a time component and is normalized: U^U^ = 1; while u M ss {7 M for non-relativistic motion of the perturbing 
matter. The perturbative nature of this ansatz is signaled by the small size of ip and X- The global behavior of the 
metric is governed by the cosmological parameters p a and p a , so that even in regions where the perturbation p is 
larger than the cosmological density p the effect of the perturbation is just to change the metric by a small amount: 

V>, X « 1. 

Einstein's equation has the form 

G„„ = -8nGT^ (16) 

where G M „ = i? M „ — (l/2)g flu R is the Einstein tensor. The spatial off-diagonal components of X^„ are small-of 0(v 2 ) 
where v << 1 is a typical non-relativistic speed of the matter in the perturbation p. With the neglect of such small 
contributions, the Einstein equation implies didjip = didjX for all i and j. This and the vanishing of ip and \ at large 
distances leads to the requirement ijj = x- The form of the time-space off-diagonal elements of G (Got — 2di(doip+-ip)) 
suggests that the substitution tp — > <p/a would be helpful, so in fact the metric that we use has the form 



^ "I 1 

The Einstein tensor for this metric, through first order in the weak field 0, is [2 




Gin/ — 
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a 2 + 2ad [ 1 - M ) _ 4ad (f> - 2aflg< 
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(18) 
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We are constrained to work to zeroth order in the flow velocity v of the perturbing mass distribution since our 
metric is not sufficiently complex to respond to the flow except through the time dependence of <fi(t,x). The metric 
perturbation <f> can be obtained from the time-time component of the Einstein equation, which is 

% ( l + -) + - ^ ^ 87rG W + p<) ■ ( 19 ) 

Here p = p m +PA represent the total cosmological energy density. We neglect the term involving the time derivative of <f> 
since d 4> is of the order of the flow velocity. Also we neglect a 2 (j> compared to V^0. Now using Mq w U 2 = g oa = 1 + 22 
and the cosmological equations ([5|), the equation for (j> reduces to 

V 2 x <j) = AnGa 3 p(t,x ) . (20) 

The "potential" 4> has the usual solution 

= - Ga J \x-x<\ = - G J (21) 

since p is the proper mass density 

p(*>* =:sr- ( 22 ) 
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Our next task is to obtain the equation of motion for a point test mass moving in the gravitational field created by 
the mass distribution p in an expanding universe. The equation of motion is the geodesic equation 

x» + r^x a aP = 0, (23) 

where the Chrostoffel coefficients are calculated from the metric of (fT7|) . They are 

r o _ F o 

1 Oi ~ 1 iO 



a \ a 

dii 

a 
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r oj = r }o = - a + -0 - d (j) I Sij , 
a \ a J 

T) k = ~ (Sikdjcj) + %9 fe 0- 5 jk di<l>) . (24) 

Then the equation of motion, in the non-relativistic and weak-source approximations x rs 1, \x\ << 1, tfi/a << 1, 
and do(j) « 1, has the form 

, < ■ d 4 . (25) 

Proceeding as for the electric case, the equation of motion in physical coordinates is 

7^--?^=-V r $. (26) 
a 

The physical-coordinate potential <j>(t, r) = 4>(t,x)/a satisfies 

= AwGp{t,r) (27) 

where p(t,r) = jJsr = ^Tpf = a? ) is the proper density of the matter perturbation in physical coordinates. For 

a stationary spherically symmetric mass distribution the potential is <& = —GM/r and (j2"B"|) has the same form as the 
electric field equation (fl"3|) . 

We note here that (|27[) . giving an instantaneous connection between mass density here and potential there, will 
break down unless the size of our local region (R) is small enough so that the cosmological mass density doesn't 
change much in the time R/c necessary for information to propagate across the region. In other words, we require 
R << Rh where Rh — c/Hq is the Hubble distance. 
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III. RADIAL EXPANSION OF ORBITS AS A MATHEMATICAL PROBLEM 

In this section we will consider the effect of the d/a cosmological force term on orbits as a purely mathematical 
problem-discussion of the physical content of this equation will be delayed until the next section. Two kinds of 
distortions of Keplerian orbits have been discussed in the literature: changes to the orbital radius and cosmologically 
induced precession. In this section we will focus on the cosmological effect on the radius. 

In order to make sense of the literature, it is essential to realize that two types of radius change have been considered. 
These are, first, the change with time of the radius of a circular orbit found for example by studying the position of 
the minimum of the effective potential; and second, the change with time of the radius of an initially circular orbit 
found by solving the equations of motion. We pause to make the distinction more clearly. 

The first type of radius change can be described more completely as follows. Let the value of the cosmological force 
term at some initial time t\ be f\. With the net force consisting of the force due to the central mass (or charge) 
plus /i, and with a given angular momentum, there is a unique radius r\ for a circular orbit. The particle will never 
actually trace out this orbit as long as the cosmological force term is changing, but this is the orbit that it would 
have were the cosmological force constant. At a later time t 2 the cosmological force has value / 2 . At this time, with 
the same angular momentum as at time f 1; there is a different radius r 2 for a circular orbit with the new force. This 



definition of the radius change Sr = r 2 — 7*1 was used by Gautreau [12j, Cooperstock, Faraoni and Vollick |13|, and by 
Bonnor [l4j in his Appendix B. The most common usage was to calculate the radius change per orbit by multiplying 
the time rate of change of Sr by the orbital period. 

The second type of radius change was found by starting a particle on what would be a circular orbit (were the 
cosmological force constant) of radius ro at time to and following the trajectory of the particle through one revolution. 
In this approach, Sr is the difference between the final radius and ro. More specifically, the particle is launched with 
the initial conditions r(to) = ro, r(io) = 0, and r(to) = 0, where ro and the angular momentum are such that, were 
the cosmological force constant at its value at time to, the orbit would be circular. The particle follows the trajectory 
required by the equation of motion, and ends up after one revolution at radius ro + Sr. This is what Bonnor [14j 
calculates in the body of his paper. Price [17j also looked at orbits found as solutions to the equation of motion. 
This solution to the equation of motion tells where the particle actually is, as opposed to the approximate position 
resulting from the first approach. 

We point out that a particle obeying the equations of motion in a time varying potential will not necessarily remain 
at the minimum of the potential energy curve if it starts there, but for a slowly varying potential the particle will 
remain at least near the minimum. This behavior will be seen for particle orbits in an expanding universe. 

In the remainder of this section we will give explicit calculations of the two types of radius change for a flat, 
dust-filled universe. We use the non-relativistic equation of motion 

r = -f + -r (28) 

r A a 

where a = GM for motion under the influence of gravitation around a central mass M; and a = —Qq/m for motion 
of a charge q with mass m around a central charge Q. For a flat dust-filled universe, the scale factor is a(t) = (t/to) 2 ^ 3 
where to is the age of the universe, so d/a — — 2/(9t 2 ). In spherical coordinates, (f2"8")) implies angular momentum 
conservation 

h = r 2 4> = const. (29) 
(where over-dots now indicate time derivatives), and the radial equation 

h 2 a 2 

^-^3+^ + ^ = 0- ( 3 °) 
A. Radius of a circular orbit 



Circular orbits have f = 0, so the radius is given by the solution to 



h 2 a 2 

7T + 3 + 7^ = 0. (31) 



For very large time the cosmological force in this model vanishes, and for a given value of h we call the corresponding 
circular-orbit radius ro. We find ro = h 2 /a. At finite times the circular-orbit radius r satisfies (|31[) and can be written 
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as r = ro(l + S). The value of r is very close to ro since the age of the universe t is so much larger than the orbital 
period r = 2irr 2 /h, which is the only other time scale in the problem. We rewrite (pJTj) as 

h 2 S--^(l + S)\ (32) 
and find that approximately S « — 2r^/ (%h 2 t 2 ). So for the orbital radius we have 



5? (i)') ■ w 



"' ' 1 ' 18, 

This is precisely result (4.10) of Cooperstock et al. [13J. It is also the result one would obtain by following the position 
of the minimum of the effective potential U(r) — h 2 /(2r 2 ) — h 2 /(r r) + r 2 /(9t 2 ) for fixed h. The effective potential 
can be written as 

™- G)'} 

where £ = t/tq. The minimum is at £ = 1 when t — > oo, but is slightly smaller than one for finite t. The position 
of the minimum grows according to p3j) as t increases. The increase of radius Q33p during the course of one period, 
evaluated at the present cosmological time to, is 

* r „* 2*(lY. (35) 
dt 9tt 2 \t J y ' 

This is equivalent to the results (6.26) of Gautreau [l2j and (B4) of Bonnor [3], but differs from (59) of Carrera and 
Guilini 0,111]. 

Of course, the result (f33|) for r(t) is not where the particle will actually be if free to move in the specified field of 
force. Instead, we have found the radius at time t of a circular orbit with a given angular momentum and with the 
cosmological force held artificially fixed. 



B. Orbital out-spiraling 



The second type of cosmologically induced radial perturbation that has been considered is that found by solving the 
equation of motion (|28p for a particle moving on an initially-circular orbit. We choose initial conditions to ensure that 
at the present moment (time to), were the cosmological force to remain fixed, the particle would move in a circular 
orbit of radius ro. That is, the initial conditions are r(to) — ro, r(to) = 0, and r(to) = 0. The time parameter t has 
an awkward zero point for our purposes here, so we reparameterize time according to 

t = t (l+s), (36) 

where the new time parameter s is small for times t not so different from the initial time to- We will use a power 
series approach to solve the equation of motion, so we write the radial variable as 



r = r 



The initial conditions immediately tell us that a\ = a,2 = 0. In fact, the vanishing of r(io) and the radial equation of 
motion (|30p give a relation among ro, h, and a: 



r 3 +3 + ^0 = 0. (38) 

' n ' n ftp 



h 2 a 2 

3 + 3 + 7)f- 
'o 'o yt, o 

In terms of the new parameters, and with use of (|3"5|) , the radial equation of motion (|3"0"|) takes the form 

d 2 n . o k — 1 2/ k 1 



, > ■ A 2 ^-^ + = ( . - .„ -^-1=0 (39) 
ds 2 k 9 V(l + s) 2 k 2 J 
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where A = ftio/rg ~ 2irto/T is a very large parameter. We find the values of the a n by equating the coefficients of the 
various powers s", and obtain the series solution for k(s): 

«M - 1 + ^ - tW + " £ V + (- H + £ V + 0( S 7 ) . (40) 



27 18 V 405 270 / V 7290 54 °, 

This series may be truncated at the s 3 term as long as s 3 << A 2 s 5 ; that is, for times t = t — to = tos small compared 
to an orbital period. For times t comparable to an orbital period, As is of order one, and the leading contributions 
for all coefficients a n with n odd are of order A™ _3 s" = 0(A -3 ). The next-to-leading terms for the odd coefficients 
have order A -5 , while the even coefficients have main contributions of order A~ 4 . This allows us to sum the principal 
contributions for t = 0(t) and get the approximate analytic form 

4 fs 3 AV A 4 s 7 \ / 1 

= 1 + JL ( As - sin(A.s)) + O (±\ . (41) 
The secular term increases linearly with time with the same coefficient 



2 



dr t 
~db ~ 9tt2 t 3 



(42) 



as in (|35p . So the actual radius follows the minimum of the effective potential as time progresses but with some 
sloshing around represented by the sin(27rf/r) term. This is in accord with expectations based on the adiabatic 
theorem of classical mechanics [29j . 

Bonnor [Til ] calculated the out-spiral radius increase starting from the relativistic equation of motion using a power 
series approach. His results (47) and thus (51) agree with our power series solution for k(s) through the term of order 
s 3 . However, in order to apply the series expansion for values of time as large as an orbital period, an infinite subset 
of terms must be summed, as shown above. We see that special relativistic kinematics do not affect the result to this 
order (as was suggested in [Til), since our calculation was purely non-relativistic. 



IV. PHYSICAL EFFECTS OF COSMIC EXPANSION 



Now we come to the real question: is there a cosmologically induced perturbation to planetary and Coulombic 
orbits, or not? In order to approach this question, we begin by considering the physical interpretation of the force. 
If the physical picture specifies the superposition of a mass M on a background homogeneous expanding universe, 
then besides the usual gravitational attraction to M there will be an additional gravitational attraction from the 
cosmological mass ^nr 3 p m inside the planet's orbit. The radial force due to this mass is f r = — G(^Trr 3 p m )^ = 
— ^f-p m r- This mass is continually decreasing on account of the outflow of cosmological mass through a spherical 
surface of radius r due to the Hubble flow, which leads to a decreasing f r . The connection of f r to the scale factor is 
through the Friedmann equations 

^r = 8^Gp, (43) 

• 2 

_^_2- = 87rGp, (44) 
or a 

where p and p are the total cosmological mass density and pressure. If we focus on a matter-only universe for the 
moment, p — p m — and ^ = -^p m . The cosmological force f r = jr in the equation of motion is present exactly 
to represent the gravitation of the cosmological mass. This agrees with the interpretation of Gautreau [T^. 

Our position then is as follows. If the model includes the cosmological mass spread homogeneously throughout the 
universe and following Hubble's expansion law, then cosmological perturbations to orbits will in fact occur. Otherwise 
there will be no perturbations to the orbits. It is imperative that on the scale set by the size of the system, it must be 
a good approximation to consider the cosmological matter to be continuous. An example of a model of this nature is 
Bonnor 's hydrogen atom situated in the Einstein-de Sitter (flat, dust-only) universe with the proton co- moving with 
the cosmological fluid [l4| . The dust must be composed of parts small compared to the size of the atom, and must 
be homogeneous as seen on a scale set by the size of the atom. In such a world, the electron (considered classically, 
of course) would spiral outwards according to the calculations of the preceding section. 
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On the other hand, if the model considered does not include actual cosmological mass expanding according to the 
Hubble law and homogeneous on a scale set by the size of the system, then there are no cosmological perturbations. 
This is the situation, for example, in the solar system. It is not possible to imagine a mass like the sun simply placed 
in a background cosmology without affecting the cosmological mass located nearby. For example, consider the sun 
and earth as a binary system. With the origin at the center of the sun, the Hubble-law speed for cosmological mass 
at the location of the earth is only 3.4 x 10~ 7 m/s, and it would take just 57fi,s for the sun to bring this mass to rest 
and initiate its plunge into the sun. Clearly a model of the solar system including cosmological Hubble flow is not 
physically reasonable. 

Now recall that the perturbation is based on the difference between the actual mass present and the cosmological 
mass. So the — V$ force in (f2"6")) is the gravitational force due to the perturbative mass. The cosmological force, 
-r, also has an associated potential — - ^- so the potential due to the total mass (cosmological plus perturbation) is 

~ -2 

$ = — + <&. It follows that the force law can be rewritten as 

Qq f , 
f=— — -V$, 45 
m r z 

where the electric term accounts for the electric force, if any. The explicit cosmological (a/a)r force term has 
disappeared: all gravitational forces are due to the mass actually present. If there is in fact cosmological mass present 
in Hubble flow, its effect is included in the — V$ force. 

The effect of dark energy is apparently a different story. There is no known dumpiness to dark energy: it seems to 
be a property of the vacuum itself, and is thus truly homogeneous. The universe at present seems to be well-described 
by a two-component model consisting of cosmological-constant dark energy and matter. Using p\ and p\ — —p\ for 
the dark energy density and pressure, and p m for the matter density (and p m = 0), the Friedmann equations become 

3a 2 

— = 8irG{p m + p A ) , (46) 
a 1 

• 2 

~ -2- = 8nGp A . (47) 
a z a 

It follows that the cosmological force is proportional to 

£ = - -irGp m + IttGpa . (48) 
a 3 3 

The p m term in a /a combines with the perturbation to give the force due to the actual, physical mass distribution. 
The dark energy term can be written as Hq£1\ where H is the present value of the Hubble constant and fl\ 

measures the density of dark energy in terms of the critical density: p\ = Qap c where p c — For values we 

use H = 70(7)(km/s)/Mpc = 2.27(0.23) x Kr^s" 1 and Q A = 0.7. ^ Including the dark energy effect (but no 
charges), the force equation takes the form 

F=-V$ + fcr, (49) 

where again $ is the potential due to the physically present matter, and kr is the linear repulsive force due to the 
cosmological constant with k = HqHa. 

The linear cosmological repulsion is a real force and is potentially observable. The challenge of course is the very 
small size of the coefficient k: k m 3.6 x 10~ 36 s -2 . The kr force is time independent, so there will be no spiraling 
of planetary orbits. The orbits will just be slightly larger than they would be without the presence of this force. 
As a modification of the 1/r 2 gravitational force law, the kr force will cause precession of elliptical orbits. We have 
calculated the rate of this precession for orbits of period r and eccentricity e, and find 

m p = ^t 2 ^t^ (50) 

47T 

radians per orbit. [3jJ This agrees with the result of Kerr, Hauck, and Mashhoon. [33| For orbits in the solar system 
this is too small to be observed. [34( An alternate expression for ([50|) is 



A0 P « 0.16 (-J-J Vl-e 2 , (51) 

which tells us that the precession is only large for systems with periods an appreciable fraction of the age of the 
universe. Excessively short human lifetime precludes any effective application of (|51|) to real orbits. 
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A potentially useful application of the cosmological repulsion is as a correction to the virial relation for galaxy 
clusters. The "moment of inertia" of a cluster can be defined to be 



I=\^M a r a -? a , (52) 

a 

where a labels the constituents of the cluster. For a cluster in steady-state, the moment of inertia is constant. This 
fact, combined with the equation of motion (|4"9")1 , yields 

= I = J2 M « (f* • K + r a ■ f a ) (53) 

a 

= M( U 2 )-^M! + 2/fc, (54) 
Ml/2 



where M is the mass of the cluster, R\i% is the radius containing half of the cluster mass, and a is a constant of order 
one defined by ^ a , b = We define an additional constant ft of order one by 21 = M (r 2 ) = (3MR\ /2 , and 



find that 

aGM 



Ri 



/2 



{v 2 )+(3kR( /2 (55) 



as the generalization of the usual virial formula. [40] As an example of the use of (|55|) . we evaluate the effect of 
the cosmological repulsion for the case of the Coma cluster discussed in Ryden Sec. 8.3 [3(|. The mean square 
speed has the (somewhat model dependent) value (v 2 ) — 2.32 x 10 12 (to/s) 2 . If we assume an exponentially falling 
mass distribution, the cosmological repulsion term contributes (3kR\^ 2 ~ (1.68)(3.6 x 10~ 36 )(4.6 x 10 22 to) 2 ks 1.3 x 

10 10 (m/s) 2 . This corrects the mass estimate by about a half of one percent, and is negligible considering the sizes 
of the other uncertainties. However, the cosmological repulsion becomes more important for more extended clusters, 
and would make a contribution comparable to that of the gravitational term for clusters that are only marginally 
bound. 

Related work includes that of Lavav et al. [Jlj], who used a cosmological constant enhanced virial relation to 
discuss galaxy formation. Current simulations of structure formation using the ACDM model include the effect of the 
cosmological constant. 



V. CONCLUSION 



Motion in local systems is governed by the gravitational forces due to the actual mass present, along with other 
forces such as electromagnetism, the linearly repulsive cosmological constant force, and initial conditions. There is 
no additional effect caused by cosmological expansion: the (a/a) f cosmological force combines with the gravitational 
force due to the perturbation of the mass distribution (mass present above and beyond the cosmological mass) to 
produce the gravitational force due to the distribution of mass actually there. It follows that calculations like those of 
Sec. 3, while interesting, do not correspond to physically realizable situations for gravitational systems. It is possible in 
principle to imagine two charged objects with one in orbit around the other, along with a continuum (at least relative 
to the size of the orbit) of out-flowing cosmological dust, for which the orbit spiraling of Sec. 3 would in fact occur, 
but there are no known systems of this type in practice. On the other hand, if "dark energy" is homogeneous with 
a constant density, the kr dark energy force is real and will affect particle orbits. We have calculated the perihelion 
precession due to the kf force, and have also found the generalized form of the virial formula relating mean-square 
speed, mass, and size for bound clusters of gravitating objects. These effects are marginally important for large 
systems such as clusters of galaxies. 
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